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ABSTRACT 

About a year after core collapse supernova, dust starts to condense in the ejecta. 
In meteorites, a fraction of C-rich presolar grains (e.g., silicon carbide (SiC) grains of 
Type-X and low density graphites) are identified as relics of these events, according 
to the anomalous isotopic abundances. Several features of these abundances remain 
unexplained and challenge the understanding of core-collapse supernovae explosions and 
nucleosynthesis. We show, for the first time, that most of the measured C-rich grain 
abundances can be accounted for in the C-rich material from explosive He burning in 
core-collapse supernovae with high shock velocities and consequent high temperatures. 
The inefficiency of the 12 C(a,7) 16 reaction relative to the rest of the a-capture chain 
at T > 3.5 x 10 8 K causes the deepest He-shell material to be carbon rich and silicon rich, 
and depleted in oxygen. The isotopic ratio predictions in part of this material, defined 
here as the C/Si zone, are in agreement with the grain data. The high-temperature 
explosive conditions that our models reach at the bottom of the He shell, can also 
be representative of the nucleosynthesis in hypernovae or in the high-temperature tail 
of a distribution of conditions in asymmetric supernovae. Finally, our predictions are 
consistent with the observation of large 44 Ca/ 40 Ca observed in the grains. This is due 
to the production of 44 Ti together with 40 Ca in the C/Si zone, and/or to the strong 
depletion of 40 Ca by neutron captures. 
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Introduction 



Massive stars (M > 8M ) are among the major contributors to the chemical enrichment 
of the galaxy through both their winds and core collapse Supernova (CCSN). In carbonaceous 
meteorites several types of C-rich presolar grains condensed in CCSN ejecta have been identified 
through association o f the measured isotopic a bundance patterns with signatures of such events 
(see e.g., IZinnerl 120031 : IClayton k Nittlen 12004 ) . Among those, C-rich grains like Type X silicon 
carbide s (SiC-X) and graph ites are assumed to require a C-rich environment to form (C/O > 



1. 



e.g., 



Travaglio et al.lll999l ). Standard massive star models predict that C-rich ejecta originate 



in explosive He shell bu rning, characterized by partial pre-explosive He-burning products (e.g. 
Wooslev k Weaverill995l ). The heavy species abundances show a mild s-process sig nature from pre- 
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explo sive convective He-shell burning (triggered by the 22 Ne(a,n) 25 Mg reaction, e.g. 
2002]), and by the explosive neutron burs t triggered again by the 2 2 Ne(a,n) 25 Mg rea ction (namely 



the n-process, with n n 



10 18 cm 3 , e.g., 



Blake k Schrammlll976l : iMever et alJ^OOOT ) 



SiC-X grains show a 28 Si excess compared to 29,30 Si, and most of these grains have a 12 C and 
15 N excess while some show a 13 C and 15 N excess. The 26 Mg and 44 Ca excess comp ared to solar 



has been explained by the condensation of unstable 26 Al a nd 44 Ti into the grains (e.g.. lAmari et al 



1992 ; iNittler k Alexander! 120021 ; iBesmehn k Hoppdl2003l ). Graphite grains (with low density, LD 



group. lAmari et al.lll995l ) show a large range of 12 C/ 13 C ratio, possibly correlated with 18 excess. 
Both excess and deficit of 28 Si with respect to 29 ' 30 Si can be observed. Indication of 26 Mg and 
44 Ca excess is also found for several of those grains. Finally, their measured N/ ° N ratio with 
few exceptions is quite normal, pointing to a contamination from solar nitrogen (e.g., IZinnerll2003l . 
and references therein). 



To explain these abundance signatures iTravaglio et al.l (I199S ) and lYoshidal (|2006l ) assume 



mixing of He-shell matter with Si/S zone material (see lMever et al 



19951 ) deeper in the supernova 



without any mixing contribution from the O-rich layers in between. The main constraint is that 
the material out of which the grains form needs to b e C-rich, to have free carbon to condense in 
SiC and graphites after CO formation. IClavtonl (|2013l ) and references therein suggests instead that 
since SN ejecta are not in thermodynamic equilibrium, C-rich grains may condense in O-rich and 
28 Si-rich ejecta. 



More recently, iMarhas et al.l ( 20081 ) found a deficit or normal 54 Fe compared to 56 Fe abundance 
in SiC-X grains , in disagreement with n ucleosynthesis predictions for the 28 Si-rich Si/S zone to be 
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54 Fe-rich (e.g ., Thielemann et al.lll996l ). Possible fractionation between Si and Fe proposed by 



(J2008j) is unlikely to be sufficient, since Fe-rich subgrains in the SiC-X grains also 



show a normal 54 Fe/ 56 Fe ratio. 

In this Letter we propose a new nucleosynthesis scenario, where the observational signatures 
of C-rich presolar grains can be obtained, without mixing, from the C-rich material in CCSN that 
experience a strong shock. The paper is organized as follows. In £j2] we describe stellar models and 
the nucleosynthesis calculations, in S|3] we compare theoretical results with measurements for SiC 
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X grains and graphite grains. Finally in $3] results are summarized. 



2. Stellar model calculations and nucleosynthesis 



In this study, we consider three SN explosion models for a 15M Q , Z = 0.02 star. (Pignatari et 
al. 2013, in preparat ion). The pre-supern ova evolution is computed with the stellar code GENEC 
(see for more details iBennett et al.l 120121. model 15S T in their Fig. 6). The explosion simulations 
include the fallback prescription bv lFrver et al.l (|2012l ). using the recommended initial shock velocity 
and reduced by a factor of 2 and 4, respectively (models 15r, 15r2 and 15r4). The sta ndard initial 
shock velocity used beyond fallback is 2 x 10 9 cms . For a 15 M Q star simulation bv lFrver et al. 
([2012]), a 2 x 10 9 cms~ 1 shock velocity corresponds to a 3 — 5 x 10 51 erg explosion, a 1 x 10 9 cms~ 1 
shock velocity corresponds to a 1 — 2 x 10 51 erg explosion, and the 5 x 10 8 cms~ 1 shock velocity 
corresponds to a weak explosion below 10 51 erg. 



The full nucleosynthesis is calculated using the post-processing code MPPNP (see e.g. . [Bennett et al 



2012 ). Here we focus only on C-rich explosive He burning ejecta which include the He/C zone and 
a small part of the O/C zone (|Meyer et al.lll995l ). 

Compared to stand ard symmetric CCSN models with explosion energy of 10 51 erg (e.g., 
Woosley &: Weaver! 1 19951 ). our models show that depending on the explosion energy 16 is depleted 
and 28 Si is accumulated at the bottom of the He/C zone and at the top of the C/O zone, where 
12 C is not significantly modified (see Fig. [2] for the details of models 15r and 15r4 with the largest 
and smallest shock). Such 28 Si-enrichment is obtained also in the intermediate model 15r2. Similar 
nucleosynthesis signatures and high tempe rature conditions ar e also predicted in high energy ejecta 
of aspherical SN and in Hypernovae (e.g., iNomoto et al.ll2009l ). 



At temperatures above ~ 3.5 x 10 8 K the (a,7)-rates of 16 and 20 Ne exceed the 12 C(a, 7)-rate 
(Fig. [I]). The reaction rate is directly correlated with the S-factor. The rate is obtained by integrat- 
ing the S-factor ove r the Gamow window, the energy range of stellar burning at a certain temper- 
ature as defined by lCaughlan et al.1 (jl985l ). The S-factor curve of 12 C(a,7) 16 is characterized by 
non-resonant reaction components: they primarily consist of interfering El transitions from tails of 
broad high energy 1~ resonances and the 1~ sub-threshold state, and a E2 dire ct capture component 
inter fering with low energy tails of 2 + resonances and the subthreshold state (jBuchmann &: Barnes 
2006). While at low temperatures the reaction rate is lower than 16 O(a,7) 20 Ne, it declines to- 
wards higher temperatures because of the decline of the S-factor by two orders of magnitude. The 
16 O(a,7) 20 Ne reaction rate on the other hand is determined by a constant non-resonant S-factor 
term and by the narrow 3~ and 1~ resonances at 1.1 Me V and 1.3 MeV respect ively, which dom- 



inate the reaction rate above temperatures of T ~ 10 9 K ( Costantini et al. 201C). Concernin g the 
reaction rates o f the subsequent a-capt ure reactions 20 Ne(a,7) 24 Mg (jSchmalbrock et al.lll983l ) and 



24 Mg(a,7) 28 Si (IStrandberg et al.ll2008l ). the non-resonant S-factor terms are larger and the critical 
energy range above 0.3 MeV is characterized by narrow resonances which cause a further enhance- 
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ment in the reaction rate, as demonstrated for 20 Ne(a,7) 24 Mg in Fig.[TJ This is also true but for 
larger energies and to a smaller extent for the production of 32 S, 36 Ar, 40 Ca, and 44 Ti. 

In parts of the He/C and C/O zones, 16 is thus depleted and feeds the production of 28 Si. 
This zone is C and Si rich, and therefore we call this zone the C/Si zone. A fundamental condition 
required to form the C/Si zone is that sufficient 4 He fuel is available while high enough temperatures 
are reached in the explosive burning. 

For model 15r, 16 and 20 Ne are depleted, whereas 28 Si and 32 S are produced up to a mass 
fraction of ~ 20% in the C/Si zone, due to reactions along the a-capture chain (Fig. [I]). As shown 
in Fig. [21 the a-capture chain is efficient up to 44 Ti. At the bottom of the He/C zone results are 
similar, but since the pre-explosive 16 O abundance is lower the production of 28 Si and 32 S is less 
efficient. In this zone Ca is depleted via neutron captures, and unstable 44 Ti is not produced. In 
model 15r4 the most produced species in the C/Si zone are 20 Ne, 24 Mg and 28 Si. In comparison 
with the pre-explosive abundances, the species above 32 S and along the a-chain are not produced 
efficiently. A similar but milder signature is obtained again at the bottom of the He/C zone. 
Notice that the rapid decrease of a-particles in the O/C zones is causing a weaker O-depletion (and 
28 Si-excess). Before the SN explosion such a zone is radiative, still carrying the signature of the 
pre-explosive convective He-burning core. The He profile under the He shell depends on uncertain 
physics assumptions in the models. Furthermore, 4 He may be mixed into the O/C zone from the 
bottom of the He shell by extra-mixing triggered by proces ses that have not been considered in 



our pre-SN model, such as rotation (e.g., iMeynet et al 



20061 ) or convective boundary mixing (e.g., 



Meakin &; Arnettl 120071 ). In that case larger fractions of the O/C zone would be converted into the 



C/Si zone compared to our model. 

In summary, in the CCSN ejecta exposed to strong shocks the upper part of the C/O zone 
and an extended region of the He/C zone are 16 0-depleted and 28 Si-rich, and 28 Si will be present 
in C-rich material when grains will start to condense. 



3. Comparison with observations 

In this section, the nucleosynthesis predictions described in ^2] are compared with observations 



for si ngle SiC-X and graphite grains, from the St. Louis Presolar Grains Database (jHvnes fc Gvngard 



20091 ). We compare only the abundances from the C-rich region in the He shell (namely, the C/Si 
zone and the He/C zone). Mildly C-rich regions in the He/N zone (H-burning layer) will be dis- 
cussed in a future work. 

In Fig. [3l Left Panel, C and N isotopic ratios for SiC-X grains are shown compared to the 
abundances in C-rich material from our models. Not included in the figure, LD graphites show a 
similar 12 C/ 13 C spread compared to SiC-X grains. The C/Si zone and the bottom of the He/C zone 
show high 12 C/ 13 C and low 14 N/ 15 N. The carbon isotopic ratio is at least 3 orders of magnitude 
larger than that of the grains, due to lack of 13 C (area labeled as "Si-rich" in the plot). Grains data 
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are located in between our predictions and the solar ratios. Some mixture with a normal component 
and /or with material from the external part of the He shell is required to fit the observations in this 
scenario. Indeed, in the most external part of the He/C zone (the part that is not O-depleted) the 
C ratio evolves to normal values, but with a N ratio larger than solar. Notice that in Fig. [3] there 
is a group of SiC-X grains showing 12 C/ 13 C lower than solar. We will discuss them in a further 
work, where SN ejecta affected by H burning are considered. In Fig.0 we report the comparison 
with LD graphites for 16 0/ 18 0, and again 12 C/ 13 C. The largest measured 18 excess is observed 
for few graphites with 12 C/ 13 C larger than solar, in agreement with the abundance signature in 
the external part of the He/C zone. In the C/Si zone and at the bottom of the He/C zone the 
16 0/ 18 becomes up to 7-8 orders of magnitude larger than solar. However, such a large ratio can 
be reduced to a ratio ~ solar with minor mixing with normal material, since there is negligible 
amount of oxygen left in these layers (see e.g., Fig. [2]). 

In Fig.|4l Right Panels, the Si isotopic ratio of stellar models and SiC-X grains are compared. 
The same is done for LD graphites in the Bottom Left Panel. In the Upper Left Panel, the 
pre-explosive and post-explosive isotopic Si abundances are reported for model 15r. During the 
explosion, in the C/O zone with low a- fuel available 28 Si is marginally produced via a-capture, and 
29,30 Si are strongly produced via neutron capture. The C/Si zone shows a strong 28 Si production. 
In particular, in zone 1 (M = 2.95 M ) the pre-explosive temperature and density are 3.1 x 10 8 K 
and 1.34 x 10 3 cm~ 3 , rising up to 2.0 x 10 9 K and 7.58 x 10 3 cm -3 respectively, duri ng the explo 



sion. Moving outward in the the deepest part of the He/C zone, the n-process (e.g., iMever et al 
2000) gradually changes the 28 Si-enrichment to strong 29 ' 30 Si-enrichments. Less extreme Si isotopic 



ratios with a mild 29 ' 30 Si-excess are seen in the external He/C zone, due to the pre-explosive s- 
processing. The 28 Si-excess is reproduced in part of the C/Si zone for models 15r and 15r2 (see 
e.g. zone 1 in model 15r, Fig.HJ Upper Left Panel and Right Bottom Panel). Also in this case, 
the grains data shows 5-values closer to the solar ratio, corresponding again to some mixing with 
unprocessed material, or with more external C-rich stellar layers. On the other hand, in model 
15r4 the shock temperature is not sufficient to reproduce the 28 Si-excess observed in SiC-X grains. 
The combined nucleosynthesis contribution from a-captures and neutron captures cause a milder 
28 Si/ 30 Si-enrichment compared to 28 Si/ 29 Si moving outward in the 28 Si-rich regions, and an earlier 
positive <5( 30 Si) with decreasing temperature compared to 5( 29 Si). This trend is consistent with a 



subset of SiC-X grains, whereas most of the SiC-X grains belong to the group XI (jLin et al.ll2002l ). 

Finally, in Fig. [J Left Bottom Panel, LD graphites show a relevant amount of grains with 28 Si 
excess, but with many scattered around the solar ratios, and few with positive Si 5 ratios. This 
different behavior between SiC-X grains and graphites is qualitatively explained looking at the 
profile of Si isotopes in model 15r (Fig. U Left Upper Panel). The efficiency of graphite condensation 
depends on the C enrichment only, whereas SiC-X grains depend on both the C and Si inventory. 
Therefore, SiC-X condensation will be more efficient in the C/Si zone than in the He/C zone, 
whereas constraints for graphites are less severe. We may also expect that the condensation of 
SiC-X and graphites will be less efficient in the external part of the He/C zone, compared to the 
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C-rich and O-depleted regions. Indeed, in these stellar layers there is less free C and Si available. 
Simulations of the condensation processes using present abundance predictions are required to 
better quantify this analysis. 

The initial 44 Ti and 26 Al en richment are distinctive signatures for a significant fraction of SiC- 
X grains and LD graphites (e.g.. IClayton &: Nittlerll2004l ), Their initial amount that condenses into 
the grains is usually estimated assuming that the 44 Ca/ 40 Ca and 26 Mg/ 24 Mg ratios are solar. This is 
not the case in our calculations (see e.g. 40 Ca abundance pr ofile in Fig.[2]) ; Furthermore, Al and Ti 
cond ense more efficiently than Mg and Ca, respectively (e.g.. lAmari et al.lll995l ; lBesmehn fc Hoppe 
20031 ). For model 15r, the thermodynamic conditions allow to produce Ca and 44 Ti at the same 
time, wit h 44 Ti/ 40 Ca ~ 10~ 2 (F ig.El). The largest measured 44 Ca/ 40 Ca ratio in SiC-X grains 
of ~ 2-3 (IHynes fc Gyngardll2009l ). would be satisfied with a fractionation Ti/Ca ~ 100 without 



considering also direct 44 Ca production. In the He/C zone, the neutron captures deplete the pre- 
explosive 40 Ca by orders of magnitude, eventually leading to high 44 Ca / 40 Ca and high 48 Ca/ 40 Ca . 
48 Ti is depleted in this case. We notice that in measurements by e.g. iBesmehn Sz Hoppd ( 20031 ) 
48 Ca and 48 Ti were not possible to distinguish. For model 15r4, the abundances in the C/Si zone 
are less extreme than model 15r, with a nucleosynthesis more similar to the He/C zone. 

One of the main puzzles for the previous ad-hoc mixing scenarios for SiC-X grains is the non- 
enhancement of 54 Fe compared to 56 Fe, since in standard CCSN models 28 Si-rich zones are also 
54 Fe-rich. In the C/Si zone of our models Fe species are destroyed by neutron captures and grain 
signatures with low 54 Fe are predicted. The isotopes 57 > 58 Fe can also be produced outward in the 
He/C zone. 



4. Conclusions and final remarks 

We have compared the abundance signature in presolar SiC-X grains and LD graphites with 
nucleosynthesis predictions from CCSN ejecta exposed to high shock velocities. Due to basic nuclear 
properties, an a-capture chain starts from 16 feeding heavier species including 28 Si, without 
affecting the pre-explosive abundance of 12 C and creating a new C/Si zone. The amount and 
distribution of a-isotopes mainly depend on the abundance of the seed 16 0, on the explosion shock 
velocity and on the amount of a-particles available. The interplay between the a-captures and the 
ra-process regulates the relative abundances of 28 > 29 > 30 Si in the He/C zone. 

The typical isotopic trends for C, N, and O of LD graphites and SiC-X grains can be reproduced 
using the C-rich ejecta, with the exception of grains with low 12 C/ 13 C. The observed 28 Si-excess 
is explained qualitatively from the signature in the C/Si zone, where 16 has been depleted by 
following a-captures to build 28 Si and other heavier species. Since production of graphites depends 
only on the C abundance, and SiC-X on both C and Si, this may explain why the 28 Si-excess is not 
observed in all LD graphites. Therefore, the C/Si zone is an ideal environment to condense SiC-X 
and LD graphites. Such abundance signatures triggered by high temperatures can be also obtained 
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in parts of asymmetric CCSN ejecta exposed to high shock velocities, or in hypernovae ejecta. 

We reconsidered the 44 Ti-excess extrapolated from a sample of C-rich presolar grains. In the 
C/Si zone the 44 Ca/ 40 Ca is not solar, affecting the results of such extrapolation. The Ca and Ti 
condensation efficiency is uncertain, leading to their elemental fractionation. For the model with 
largest explosion energy (model 15r) both 40 Ca and 44 Ti (and 48 Ti) are produced in the C/Si zone 
(with 40 Ca/ 44 Ti ~ 100). Neutron captures destroy Ca in the He/C zone. The isotope 44 Ca 
(and 48 Ca) instead is produced by neutron captures in all the O-depleted and C-rich ejecta. For 
weaker shocks (see e.g., model 15r4) neutron captures dominate the Ca nucleosynthesis, 40 Ca is 
depleted and 44 Ca produced. In both cases, a compari son needs to be made gr ain- by- grain, where 
correlations with e.g., Ti species and 28 Si-excess (e.g.. IClavton fc Nittlerl 12004 ) are considered. A 
similar discussion can be derived for 26 Al enrichment. 

Finally, the present results provide support to the existence of different energy components 
in the ejecta of the same CCSN. Asymmetries during the SN explosion of e.g., SN1987A and 
CasA are confirmed by observations. In SN1987A the emission Sil and Fel-II from core material 
are correlated with Hel emission, possibly indicating the signature of q-rich freezo ut during high 
energy -SN explosion in at least part of the ejecta (jKiaer et al.ll2010l ). In CasA, iDeLanev et al 



(J201Cj) reported two major velocity components for Sill in the ejecta. We suggest that they may be 
the signature of the Si/S-Si/O and C/Si zones, respecti vely. This scenario would also explain the 
non-correlation between Nell and OIV emission in CasA (jlsensee et al.ll2010i ). Indeed, Ne (together 
with Mg, S, Ar and Ca) could be produced also in the C/Si zone, whereas for standard CCSN O 
and Ne are produced in the same regions. 
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Fig. 1.— The S-factors for the 12 C(a,7) 16 0, 16 O(a,7) 20 Ne and 20 Ne(a,7) 24 Mg are shown as 
function of a-energy in the laboratory system (Left Panel). The a-capture rates of the a-chain 
from 16 to 44 Ti and normalized to the rate of 12 C(a,7) 16 are shown as a function of temperature 
(Right Panel). 
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Fig. 2. — Abundances profiles in the C-rich material for H, 4 He, 12 C, 16 and 28 Si before (thin 
lines) and after (thick lines) the supernova explosion for models 15r and 15r4 (Upper Left and 
Right Panels, respectively). The same is shown in more detail and for more species in the lower 
panels. The C/Si zone is where O is depleted producing 28 Si, above 2.945M©. The process of 
O-depletion affects a large part of the He/C zone, up to M = 3.3M and 3.2M Q for models 15r 
and 15r4, respectively. 
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12C/13C 12C/13C 

Fig. 3. — Isotopic ratios 14 N/ 15 N and 12 C/ 13 C from models 15r, 15r2 and 15r4 compared with 
measured data for SiC-X grains (Left Panel), and 16 0/ 18 for LD graphites (Right Panel). The 
numbers in the plots are mass coordinates where the C-rich region starts and where it ends in the 
ejecta. Different points along the lines indicate different mass-zones where the nucleosynthesis was 
calculated. Continuous black lines indicate the solar ratios. 
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Fig. 4. — Si abundances in the Si-rich material for model 15r before and after CCSN (Upper Left 
Panel). Si isotopic ratios in 5 notation (J(ratio) = (stellar ratio/solar ratio - l)xl000) is given 
for SiC-X grains (Upper Right Panel and zoomed in the Lower Right Panel) and LD graphites 
(Lower Left Panel) compared with models. The Si abundances in model 15r and the correspondent 
5-values are highlighted for 4 zones, at M = 2.95, 3.15 (zones 1 and 2), 3.3 and 4 M Q (zones 3 and 
4, out of the x-axis range in the Upper Left Panel). 



